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A phenomenological thermodynamic theory of Ba^iSri-^TiOa (BST-a;) thin films epitaxially 
grown on cubic substrates is developed using the Landau-Devonshire approach. The eighth-order 
thermodynamic potential for BT single crystal and modified fourth-order potential for ST single 
crystal were used as starting potentials for the end-members of the solid solution with the aim to 
develop potential of BST-x solid solution valid at high temperatures. Several coefficients of these 
potentials for BT were changed to obtain reasonable agreement between theory and experimental 
phase diagram for BST-x [x > 0.2) solid solutions. For low Ba content we constructed the specific 
phase diagram where five phases converge at the multiphase point {Tm2 = 48 K, a; = 0.029) and 
all transitions are of the second order. The "concentration- misfit strain" phase diagrams for BST-a; 
thin films at room temperature and "temperature- misfit strain" phase diagrams for particular con- 
centrations are constructed and discussed. Near Tn2 coupling between polarization and structural 
order parameter in the epitaxial film is modified considerably and large number of new phases not 
present in the bulk materials appear on the phase diagram. 

PACS numbers: 64.60.Kw, 64.70.Kb, 77.65. Ly, 77.80.Bh, 77.84.Dy, 81.30.Dz 



I. INTRODUCTION 

Ferroelectric thin films are very good candidates for a 
wide range of applications: for example, as high-density 
dynamic random access memories, large-scale integrated 
capacitors, pyroelectric detectors, micro- and nanoelec- 
tromechanic devices Epitaxially grown thin films 

are usually highly constrained because their fabrication is 
accompanied by several strain factors originating due to 
lattice mismatch and the difference between the thermal- 
expansion coefficients of the film and the substrate. Also, 
transformation strain usually appears at the ferroelectric 
phase transition if the heterostructure deposed at ele- 
vated temperatures is cooled below the Curie tempera- 
ture (Tc). The importance of strain effects in determining 
the properties of thin-film ferroelectrics was recognized 
in numerous theoretical and experimental investigations 

siiiBii[ii[ii[ii,[a. 

The Baa;Sri_a;Ti03 (BST-a;) solid solution is one of 
the most studied lead-free ferroelectric system from fun- 
damental and technological point of view. In the para- 
electric cubic phase, the end-members BaTiOa (BT) and 
SrTiOa (ST) have similar crystal structure but exhibit 
quite different phase transition sequences on cooling. 
BT displays three phase transitions towards ferroelectric 
phases characterized by the polarization along different 
directions. On the other side, ST shows no ferroelectric 
phase transition on cooling but first an antiferrodistor- 
sive transition followed by the appearance, at very low 
temperature, of a quantum paraelectric state. Both BT 



and ST compounds have been widely investigated, both 
theoretically and experimentally, whereas, the efforts on 
the BST-x solid solution are mainly experimental. It is 
then very desirable to develop new calculation and the- 
oretical tools that allow better and deeply describe the 
structure and thus the properties of this system. The 
temperature-concentration phase diagram of bulk BST- 
X solid solutions based on the phenomenological six-order 
thermodynamic potential was recently developed [141] and 
compared with available experimental data. 

Thin films of BST-x solid solutions are of great tech- 
nological interest due to their excellent ferroelectric and 
piezoelectric properties. The structure and properties of 
ferroelectric thin films substantially differ from those of 
bulk ferroelectrics; therefore, theoretical analysis is im- 
portant since it can provide fundamental insights into 
the behavior of thin films. Recently, Ban and Alpay [lH| 
have developed phase diagrams for single-domain epitax- 
ial Bao.6Sro.4Ti03 and Bao.rSro.sTiOa films on (001) cu- 
bic substrates as a function of the misfit strain based on 
the Landau-Devonshire six-order potential as a polyno- 
mial of the polarization components in accordance with 
the one developed by Pertsev et al 0]. The parameters 
used for the calculation of the renormalized coefficients 
for BST films were obtained by simple averaging the cor- 
responding parameters of BT and ST. It is worth noting 
that the contribution of six-order polarization terms to 
the free energy was neglected in their calculations that is 
rather oversimplified approach. However, as known from 
experimental investigations, high-order terms should be 
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taken into account, because the two-dimensional clamp- 
ing increases considerably the ferroelectric phase transi- 
tion temperature in perovskite thin films. Stress-induced 
upward shift in Tc as large as tens of degrees have been 
observed In this case, the coefficients at the sixth- 

order terms in the thermodynamic potential 5] becomes 
negative, therefore, the six-order expansion is not valid 
for high-temperature phase transitions in thin films. Very 
recently , the eighth-order Landau-Devonshire poten- 
tial was used to construct " misfit-temperature" phase di- 
agrams of epitaxial BT thin films on cubic substrates. 

The present paper is devoted to BST-x thin films and 
is organized as follows. First, we develop thermodynamic 
potential for BST-a:; solid solutions using phenomenolog- 
ical models known for pure BT and ST single crystals. 
Due to the above mentioned reason the sixth-order po- 
tential for BST solid solutions used in [14] is not valid 
for thin films, therefore here we use eighth-order po- 
tential developed by Li [13] for BT single crystal and 
fourth-order potential recently developed for ST crystal 
p^ . It is worth noting that in [13] the coefficient at 
in the potential for ST was changed considerably to 
achieve agreement between experimental data 19] and 
para-ferroelectric transition line on the theoretical phase 
diagram. As a result, the value of the susceptibility was 
found to be significantly and abnormally large with re- 
spect to the experimentally observed one. In the present 
work we use coefficients of thermodynamic potential for 
ST [11] and BT [l^ Moreover, several coefficients (Qn, 
Qi2j 0:123) from [i3] were changed to get better agree- 
ment with experimental diagram in the Ba-rich side. Sec- 
ond, we develop thermodynamic potential for BST thin 
film epitaxially grown on cubic substrate. The potential 
derived for BST solid solutions according to the method 
proposed in [31 includes no term associated with ther- 
mal expansion. To take into account thermal expansion 
the above method was modified with the aim to include 
evidently relevant term in the potential. Finally, the re- 
sulting potential was used to construct according to the 
conventional way the "misfit-temperature" diagrams 
for thin films of particular BST-x compositions. 



II. A PHENOMENOLOGICAL 
THERMODYNAMIC POTENTIAL FOR BST 
SOLID SOLUTION 

Following the method developed in ^li"] the Hclmholtz 
thermodynamic potential of the solid solution can be 
written using the known thermodynamic potentials 
FsriVi^) ^iid FBTiVi^) for end members of the solid 
solution x = (ST) and a; = 1 (BT): 

F = (1 - x)Fst{v, u - Ast) + xFbt{v, u - Abt) (1) 

where 77 is the order parameter, m is the common elas- 
tic strain of the solid solution, Ast and Abt are lattice 
strains of the end members, needed to fit the lattice pa- 
rameter ST, BT and BST- a; Ti?]. The lattice parameters 



of the solid solution end-members are 



asT - 
a-BT 



Ast), 
' Abt)- 



(2) 



The lattice parameter of the solid solution ax can be 
found from the condition, which implies complete com- 
pensation of internal elastic forces: 



(1 



du 



ri,u—0 



ri,u—0 







(3) 

The thermodynamic description may be developed 
starting from the power-series expansion of the Gibbs po- 
tential $(ry,i). The relevant general expression for cubic 
perovskite is given as 



+aii {p\+pi+pt) + "12 {pIpI + pIpI + pIpI) 

-hi {flpl + flpl + fhl) - ti2[fl{pl +pI) 

+ ipl{pl +pI) + (pUpI +pI)] ~ Uiiip2'P3P2P3 
+ ipiip3PlP3 + V1V2P1P2) -I- $6 + *8 + ^'t 

(4) 

where p is the order parameter - polarization related to 
ionic shifts in polar zone-center mode, Lp is the out- 
of-phase rotation of TiOg octahedra corresponding to the 
R25 zone-boundary mode in the cubic phase Pm'im(0\) . 

The high-order terms are written as 



pI) 



pUpI 



pI) 



*6 = am {pt + p^ + pf) + aii2[ptipl 

+Pt{pl +pI)] + ai23Plplph 
$8 = aiiii(pf +Pi +P3) + aiii2[p?(pi + pI) + pUpI + pD 

+pI{pI +pI)] + aiMpjpt+ptpt+ptpi) 

+aii23{p\plpl + plpipl + pIpIpI) 

(5) 

The elastic energy $t is: 



*t = -^11 {ti(p\ + t2(p\ + hvl) - Rn{tAV2V3 + hVW3 
+t6^1^2) - Rl2[tl{^l + (fit) + t2{Lpl + ifil) 

+ ipl)] - On [tipl + t2pl + hpf) 

-Qn lUP2P3 + tbPlP3 + hPlP2) 

-Q12 itlipl+pl)+t2ipl +Pl)+t3ipl +pI)) 

-5SII (tl +4+ tj) - is44 (tl +tl+ tl) 

-si2 {ht2 + ht-3 + t2h) - XT {tl + <2 + h) 

(6) 

where Skj are the compliances, A is the linear thermal 
expansion coefficient, T is the absolute temperature, the 
stresses ti, i=1..6 are given in Voigt notations. In the 
following, the coefficients of the potentials (4)- (6) of the 
end-members of the solid solution are denoted by addi- 
tional indexes ST and BT. The reference point on the 
temperature scale is taken at T = 0. Temperature T 
should be substituted in Eq. 6 by AT = T - To if refer- 
ence point To is not equal to zero. 

Eq. (1) is written for Helmholtz potential F{r},u) , 
which can be obtained from Gibbs potential (4), (5), (6) 
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by the formal substitution of the coefficients and sub- 
stitution of stress t in (4) for the quantities determined 
from the equations Ui = [lO]- Only coefficients at 

second and fourth order terms in (4) will be renormalized 
because Eq. (6) contains quadratic strains. 

The Helmholtz potential of the solid solution (1) con- 
tains no distinct term corresponding to the thermal ex- 
pansion. Actually, temperature dependence of the lat- 
tice parameters is included in the lattice parameter of 
the solid solution Ux- To construct the potential in a 
conventional form, when linear thermal expansion is as- 
cribed by the terms linear with respect to strains, one has 
to shift the common strain u by the value of the linear 
thermal expansion. The latter can be found from the lin- 
ear expansion of (Eq.7 in Ref.l4). However we used 
the following way. In Eq.3 we exclude the terms cor- 
responding to thermal expansion, therefore we assume 
AgT and Agx to be temperature independent. In this 
case, AsT and Abt derived from Eqs. (2)- (3) preserve 
the term corresponding to the thermal expansion in Eq. 
(1). Under the above introduced conditions, in the cubic 
paraelectric phase Eqs. (2) and (3) yield: 

— x6 



AsT — 
Abt — 



(1- 



+ (l+<5)a 
x\r5 



(l-a:)T + (l-|-i5)a; ' 
{l—x)TasT+xaBT 
{1 — x)t+x 



(7) 



The thermal expansion coefficient can be written as: 

((1-.x)t+(1 + 7)x) ^ 
A.T — ^ AsT \o) 



{1 — x)t + X 



where r 



Eq. (8) can be easily obtained from Eq.7 of Ref. 14 if 
the lattice parameter of the solid solution is presented as 
one-series expansion ubst ^ CLxi^ + ^xT). 



' asT 



III. PHASE DIAGRAM OF THE BST-x SOLID 
SOLUTION 

The starting Gibbs potential for pure BT is expanded 
as pol yno mial of polarization components up to eight 
order [17| and contains only additional ip"^ term with 
the positive coefficient equal to 3.7 x 10^^ J/m^ Ref. 
14. All coefficients are listed in Table I, where coeffi- 
cients in Eqs. (4)- (6) in higher-order terms containing 
if are equal to zero. Note, the coefficient ai23 used in 
Ref. 17 was changed in our calculations to fit theoreti- 
cal orthorhombic-rhombohedral phase transition temper- 
ature to the experimental value. All coefficients in Gibbs 
potentials for end-members BT and ST used in this work 
are listed in Table I. 

It is important to emphasize that it is not easy to 
deal with the coefficients especially in the Sr-rich re- 
gion where both the quantum effects and the coupling 
between ferroelectricity and antiferrodistorsivity are co- 
existing. As a consequence, the coefficients of the poten- 
tial depend strongly on the technological parameters of 



TABLE L Coefficients of Gibbs potentials for BT [17| and ST 
[l^ . Coefficients in the quadratic terms are: /3i,st ~ 1.036 x 
10^®[coth(43.8/r) - coth(43.8/106)] , /3i,bt = 3.7 x 10^^ in 
J/m^ Qi,sT = 4.05 X lO'^ X [coth(54/r x ) ~ coth{54/30)], 
ai,BT = 4.124 X 10^(r- 388) in Jm/C^ Coefficients from 
Refs. [itI . [21I ] different from those used here are given in 
brackets. 
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xlO-2 mVC^ 


Q12 


-1.31 


-4.5 




Q44 


1.9 


2.9 




Sll 


3.52 


8.33 


xlO"^^ mVJ 


S12 


-0.85 


-2.68 




S44 


7.87 


9.24 





the samples synthesis including homogeneity, impurities, 
vacancies, stress conditions, grain boundaries, but also 
on the techniques (Raman spectroscopy, dielectric spec- 
troscopy, x-ray diffraction, etc) used to have access to 
these coefficients. As known, even dielectric properties 
are very sensitive to the conditions of the sample prepa- 



ration 
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thus coefficients in the terms depend on 
these parameters too. Also, as follows from the proposed 
method of derivation of thermodynamic potential of solid 
solution, any impurities induce additional deformations, 
which renormalize coefficients of the potential. As a con- 
sequence, coordinates of critical points can be different 
for samples in single crystal, powder or ceramic form 
prepared in different laboratories. Nevertheless, even if 
the critical points may differ from one sample to another 
one, the qualitative phase diagram should not change and 
can serve as a starting setting to study the consequences 
when BST-x is as thin film form. 

The resulting phase diagram for BST-x solid solutions 
calculated using the coefficients listed in Table I is pre- 
sented in Fig. 1. Although phase diagram has very simi- 
lar overall view with respect to that previously reported 
[l^ . there are some changes caused by the above men- 
tioned changes of the coefficients in fourth-order terms: 
1). The A^-phase point T^i is now shifted from x = 0.13 
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FIG. 1: Phase diagram of bulk BST-s solid solutions. Solid 
lines and dashed lines correspond to first and second order 
phase transitions, respectively. Experimental points (circles) 
taken from Ref. 19. The calculated coordinates of the tri- 
critical and multicritical points are Tjvi = 116 K,2; = 0.2; 
Tn2 = 48 K,x = 0.029; = 317 K,a: = 0.75. Detailed 
diagram at low x is presented in Fig. 2 
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FIG. 2: Calculated phase diagram of bulk BST-s in the vicin- 
ity of multiphase point Tn2 = 48 K, a:: = 0.029. Dashed lines 
correspond to second order phase transition 



to a; = 0.11. 2). The tricritical point now appears at 
X = 0.37. 3). At low X five possible diagrams were 
proposed in |14l | for different possible i44 values. Since 
coefficient was determined [18] from the experimental 
data, only one diagram at low x values is calculated (see 
in Fig. 2). 



IV. PHASE DIAGRAMS OF THIN FILMS 

Below we consider BST-x thin films epitaxially grown 
onto (001) surface of a cubic substrate. Deposition usu- 
ally occurs well above the of bulk BST-x. The film 
is constrained due to lattice mismatch between the film 
and the substrate. This deformation occurs in the plane 
parallel to the substrate and induces symmetry lower- 
ing from cubic Pm3m to tetragonal PA/mmm. As 
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FIG. 3: The "concentration-misfit strain" diagram for BST-a; 
thin films at T = 3QQK. The coordinates of the multiphase 
points are Ti(u™, = 0, a: = 0.7), 72 (Wm = -15.5 x 10^, a; = 
0.083) and Ti{uru = 0.26 x 10^, a; = 0.81). Dashed lines corre- 
spond to the second-order phase transitions. Phase number- 
ing is given in Table II. 



a consequence, all three-dimensional irreducible repre- 
sentations of the cubic phase split into one- and two- 
dimensional representations. The group-theoretical anal- 
ysis of the structures derived from the aristotype cubic 
perovskite Pm3m by two order parameters (f i„ and R25 
soft modes) yields 26 low-symmetry phases, while sym- 
metry lowering of the paraelectric to tetragonal phase in- 
duces 33 phases. A list of possible low-symmetry phases 
derived though TiOe octahedral tilting and Ti-cation dis- 
placements is given in Table II together with their nota- 
tions already used in previous literature [l^l ■ 

In the absence of external forces the strains at the 
film substrate interface are un ~ U22 ~ Um, U12 — 0, 
where the misfit strain u,„ exhibits temperature depen- 
dence [i6| : 

as 



7.0 + (A - -^Xs) (T - To) 



where uo 



S11+S12 a—c 



is the primary deformation of the 



film at the deposition temperature Tq, a and c - the cor- 
responding in-plane and out-of-plane lattice parameters 
of the film at Tq, 05 is the lattice parameter of the sub- 
strate, A and As - the thermal expansion coefficients of 
the film and the substrate, respectively. 

Following Pertsev et al [5[ we use the Gibbs thermody- 
namic potential (4) with the coefficients listed in Table I 
to construct the " concentration- misfit strain" [x — u„i ) 
phase diagrams for BST thin films at room temperature 
and "temperature-misfit strain" (T — Um ) ph ase diagrams 
for particular concentrations. Recently [la. [25j. we have 
illustrated that phase diagrams for epitaxial BT films 
depend on the values of compliances and electrostrictive 
coefficients used in calculations. In the present work we 
performed calculations using the set of coefficients suit- 
able for the whole phase diagram of BST-x solid solu- 
tions. 

The X — Um phase diagram for BST-x thin films at 




FIG. 4: Um — T phase diagrams of BST-a; thin films for selected concentrations. Solid lines and dashed lines correspond to 
first and second order phase transitions, respectively. Phase numbering is given in Table II. Detailed diagrams at low x around 
Um = is presented in Fig. 5. 



300 K is shown in Fig. 3. There are five low-symmetry 
phases and two multiphase points Ti; (u„i = 0, x = 0.72), 
Ta: (li™ = -15.7 x 10^, a; ^ 0.1) and T3: (u„ = 
0.26 X 10"^, a; = 0.81), which is very sensitive to the Qij 
values. The tetragonal paraelectric phase may be stable 
at room temperature up to x = 0.72 at zero misfit strain. 
The monoclinic phase 5 - {pi ^ 0, p2 = and ^ 0) , 
with the polarization tilted with respect to the film sur- 
face is stable above Ta point and at low misfit strains 



The tetragonal phase 1 with the polarization normal to 
the substrate {pi = P2 = and p^ ^ 0) is stable in a 
wide range of negative misfit strains and x > 0.1. In 
the right side of the diagram the orthorhombic phase 3 
( Pi = P2 7^ and Pa = 0) is stable. These phases were 
already discussed in previous literature [l^ for BST-x 
thin films. As follows from Fig. 3, new phases may ap- 
pear at room temperature only at high negative misfit 
strains and only in the Sr-rich side of the phase diagram, 
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TABLE II: Low-symmetry phases allowed for BST-x thin 
film deposited on (001) cubic substrate. For all phases with 
nonzero order parameter ip the primitive cell volume is dou- 
bled as compareed to that in the high-symmetry phase. The 
i/Ji and Pi components correspond to _R25 and F\u soft mode, 
respectively 
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where the two-dimensional (2D) clamping (negative Um) 
stabilizes the I A/ mem tetragonal phase 8, which exists 
in bulk ST crystals below 106 K. Finally, very strong 2D 
clamping stabilizes tetragonal lAcm polar phase 16 in 
the narrow concentration interval in the vicinity x ~ 0.1. 

Three-dimensional u„j ~ x~T phase diagram is rather 
complicated and contains the multiphase point T/v (um = 
0,x = 0.029, T = 48 K where all second-order terms in 
Gibbs potential are equal to zero and a large number of 
low-symmetry phases converge. Below we discuss several 
Um — T diagrams for particular concentrations. Phase di- 
agrams for X = 1.0 and 0.8 presented in Fig. 4 are similar 
to one of the diagrams reported for BT in Ref. 16. In 
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FIG. 5: Phase diagrams of EST- a: thin films at low x around 
Um ~ 0. Solid lines and dashed lines correspond to first and 
second order phase transitions, respectively. Phase number- 
ing is given in Table II. 



the narrow interval between paraelectric phase and or- 
thorhombic phase 3 with in-plane polarization along the 
basal diagonal of the unit cell there is phase 2 with the 
polarization along one of the former cubic axis in plane 
of the film. Also, in the narrow interval between phase 1 
with the polarization normal to the substrate and phase 
4 with polarization along the space diagonal of the unit 
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cell there is monoclinic phase 5. The stability ranges of 
these intermediate phases (2 and 5) decreases with the 
decreasing of Ba content and below x — 0.74 they are not 
stable any more. As a result, phase transition 0^3 and 
1 <-> 4 are the second order, however, there is rather nar- 
row temperature interval where first order 1 <-> 3 phase 
transition occurs (see diagram for x = 0.6 in Fig. 4). Be- 
low X = 0.4 this first-order boundary disappears and only 
second-order phase transitions occur for x = 0.4 and 0.2. 
In this concentration interval phase diagram is similar to 
that developed for the fourth-order potential. Only three 
low-symmetry phases (1, 3 and 4) present on the phase 
diagrams and converge in the multiphase point Ti , which 
steadily decreases with increasing Sr content and obeys 
the same law as the paraelectric-tetragonal phase tran- 
sition line in Fig. 1. For x = 0.08 one more multiphase 
point T2 appears on the left side phase diagram and two 
additional phases, namely 8 and 16, are stable for nega- 
tive misfit strains. Tetragonal phase 8 {QOcp 000) is sim- 
ilar to that in bulk ST crystal below 106 K, while phase 
16 {OOip OOp) is induced by the mixed order parameter 
including both octahedral tilting and ionic displacement 
in the direction normal to the substrate. 

The multiphase point T2 moves towards Ti with de- 
creasing X and these points coincide at x = 0.029. Be- 
low this critical concentration, the T2 point appears now 
on the other side of the diagram at Um > 0, while Ti 
point is always at Um = 0. Three diagrams for x=0.025, 
0.01, and are presented in Fig. 4. Two low-symmetry 
phases 8 {OOip 000) and 9 {ipOO 000) converging at Ti 
correspond to different domains of bulk ST. The transi- 
tion line between these phases is of first order. In Fig. 
4, the paraelectric phase and five low-symmetry phases 
converge in the T2 point. According to Ref. 14, the 
existence of such multiphase point r2, where five low- 
symmetry phases converge, does not contradict the Gibbs 
phase rule. Below critical concentration x = 0.029, both 
Ti and T2 move to higher temperatures with decreasing x 
and below T2 phase diagrams in Fig. 4 change consider- 
ably. Namely, six phases 0, 3, 15, 32, 21 and 9 converge 
in the multiphase point T2. 

Five phases allowed below 60 K for a; = are shown in 
Fig. 5c. Phase transitions 8^9 and 18 ^21 are of the 
first order, while 8 ^ 18, 9 ^ 21 and 21 ^ 32 are of 
the second order. Note all these phases do not converge 
in one point. There are two close tricritical points so that 
phase transition line between 9 and 18 exists. New phase 
27 appears on the phase diagram at x =0.01 in between 
18 and 8 phases, while phase 16 is stable for negative 
misfit strains (see Fig. 5b). 

Near the critical concentration x = 0.029 the phase 
diagrams are very complicated. Coupling between polar- 
ization and structural order parameter in the epitaxial 
film is modified considerably and new phases that were 
not present in the bulk material (Fig. 2) appear. Fig. 4 
shows overall view and Fig. 5a shows detailed diagram 
around Um = for x = 0.025. There are thirteen low- 
symmetry phases in a film instead of six in a bulk sam- 



ple. Both rhombohedral phases i?3c and i?3m allowed in 
a bulk material are not stable in the film. 

Finally, we compare our phase diagram for pure _ST 
film (a;=0) with those developed by Pertsev et al |24|. 
Although phase diagrams have similar overall view for 
Um X 10^ |< 10, some important differences caused by 
the above discussed changes of the coefficients should be 
emphasized. As shown in Fig. 4, additional phase 32 ex- 
ists in between 21(F0II) and 15(F0IV) phases. We have 
to emphasize that FOIV phase determined as | pi | = | P2 | 
and I V'l I = I '/'2 I (Ref. 24) actually contains two different 
phases 15 and 17 listed in Table II. Note, that ferroelec- 
tric (FTI) tetragonal phase 1 (000 OOp) predicted in Ref. 
27 for negative misfit strains and T >200 K, according to 
our calculations is stable at low temperatures only (see 
Fig. 4). Also, phase 15(F0IV) exists in the limited misfit 
interval. 

It is worth noting that in pure ST many phases are sta- 
ble in the narrow intervals near Wm= 0. Therefore, even 
weak clamping of single-crystalline plates, usually used 
for low-temperature measurements, can induce phase 
transition to one of the low-symmetry phases. 

Concluding this section we have to emphasize some 
limitations of the above developed approach. The phe- 
nomenological potential of the solid solution (1) was writ- 
ten through the potentials of its end members {x = 0, 1). 
Therefore the validity of the potential (1) is based on 
validity of the potentials for each member of the solid so- 
lution. The method used to develop the potential (1), (2) 
and (3) is based on the macroscopic theory of elasticity 
and crystal lattice matching. We assume that solid solu- 
tion with no ordering in the whole concentration range 
can be correctly described if ionic substitution induces 
geometrical distortions, accompanied by the macroscopic 
elastic deformation [26| . We also assume that isomorphic 
Ba/Sr substitution does not change the type of chemical 
bonding. Therefore, in our model, most (or even all) 
changes in the solid solution are caused by the elastic 
interaction. 

Actually, currently available potentials for both end- 
members of BST solid solution require further improve- 
ments. Some coefficients for fourth-order potential devel- 
oped by Uwe [13] for ST were recently revised [ll|. Al- 
though some uncertainty in the values of the coefficients 
still exists, reasonable agreement with the available ex- 
perimental data were found ^18'|. Further investigations 
and new perfect experimental data will require extension 
of the theory to the higher orders [2^ . 

Two eight-order potentials are currently available for 
BT. Li et al [l7|] developed a model with only a single 
temperature-dependent second-order coefficient, while all 
higher-order coefficients were temperature-independent. 
This model could successfully reproduce the phase transi- 
tion temperatures and their dependence on electric fields, 
as well as the dielectric and piezoelectric constants of 
bulk BT. Recently, Wang et al [1^ developed poten- 
tial where forth-order coefficients were also temperature- 
dependent and successfully explained dielectric proper- 
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ties of BT. Although the latter potential looks much more 
attractive, we failed to employ it in our model for BST 
solid solution. Since elastic properties of the potential are 
very important in our model, we have analyzed the tem- 
perature - hydrostatic pressure (T — p ) phase diagram 
using potential [2§|. We found that the paraelectric- 
ferroelectric phase transition temperature decreases with 
increasing pressure, approaches the minimum, then in- 
creases, and finally decreases again. Such a behavior 
contradicts to the experiment [30|. Although the po- 
tential developed by Li does not reproduces exactly 
the experimental T — p phase diagram, a disagreement is 
not so pronounced. Very likely, as was noted in Ref. 31, 
additional terms should be considered in the potential to 
describe large deformations. 

V. CONCLUSIONS 

To construct the " concentration- misfit strain" and 
"temperature- misfit strain" phase diagrams of epitaxial 
BST- a; thin films on (OOl)-oriented cubic substrates we 
modified thermodynamic potential for bulk BST-x solid 
solutions developed in Ref. 14. The thermal ex pan sion, 
missing in the previous thermodynamic model |14| was 
simply introduced in the Helmholtz potential (1) by the 
shift of the common strain u of the solid solution on the 
value of the linear thermal expansion. The introduced 
modification docs not change the phase diagram of bulk 
BST- a; solid solution, and allows ones to develop thermo- 
dynamic theory for thin films. 

In the present work we used eighth-order thermody- 
namic potential for BT single crvst al [l7| and fourth- 
order potential for ST single crystal [I^ to develop rel- 
evant potential of BST-x solid solution, which can be 
applied at high temperatures. 

In order to fit the phase transition lines to experimen- 
tally measured phase diagram of BST-x solid solutions 
p^ . the coefficients Qii and Q12 in six-order term in the 



potential of BT were changed. Also, the stable thermody- 
namic states and phase transitions lines in the concentra- 
tion range x < 0.04 were calculated. Five low-symmetry 
phases are predicted for bulk solid solution in the narrow 
concentration interval around multiphase point Tjv2 = 
48 K, X = 0.029 as shown in Fig. 2. 

Performing the necessary calculations, we constructed 
the room-temperature " concentration- misfit strain" {x — 
Um) phase diagram for BST-x thin films epitaxially 
grown on cubic substrate. The diagram is useful for prac- 
tical applications in thin-film engineering. Depending on 
the type of strain imposed by the substrate, epitaxial 
BST- a; thin films can be grown with the polarization 
either normal or parallel to the substrate. As follows 
from the phase diagram presented in Fig. 3, ferroelectric 
state with the polarization parallel to the substrate is 
only possible is ST films deposited on tensile substrates. 
Recently [32| , room-temperature ferroelectricity was ob- 
served in ST films on DyScOs substrates. The tensile 
lattice mismatch in this heterostructure at room temper- 
ature is about 1% that is very close to the 0-3 phase 
transition line in Fig. 3. 

A rich variety of low-symmetry phases induced by oc- 
tahedral tilting and ionic displacements are allowed near 
the TV-phase point (m™ = 0,x = 0.029, T = 48 K). Most 
of them are stable in the narrow temperature intervals 
near Um= 0. Great difference between the phase tran- 
sition sequence in bulk BST-a: and epitaxial thin films 
with the composition around x — 0.029 can be therefore 
expected. The selection of the substrate and the film 
composition allow manipulating the strain state in the 
film to achieve the desirable phase transition sequence. 
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